Bolometric technique for high-resolution broadband microwave spectroscopy of 

ultra- low- loss samples 
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A novel low temperature bolometric method has been devised and implemented for high-precision 
measurements of the microwave surface resistance of small single-crystal platelet samples having very 
low absorption, as a continuous function of frequency. The key to the success of this non-resonant 
method is the in-situ use of a normal metal reference sample that calibrates the absolute rf field 
strength. The sample temperature can be controlled independently of the 1.2 K liquid helium 
bath, allowing for measurements of the temperature evolution of the absorption. However, the 
instrument's sensitivity decreases at higher temperatures, placing a limit on the useful temperature 
range. Using this method, the minimum detectable power at 1.3 K is 1.5 pW, corresponding to a 
surface resistance sensitivity of «1 /i£7 for a typical 1 mmxl mm platelet sample. 



I. INTRODUCTION 

In this article we describe an apparatus designed for 
the continuous-frequency measurement of low tempera- 
ture electromagnetic absorption spectra in the microwave 
range. The motivation to develop this instrument comes 
from a desire to resolve, in great detail, the microwave 
conductivity of high-quality single crystals of high-Tc 
cuprate superconductors. However, the technique we de- 
scribe should find a wealth of applications to other con- 
densed matter systems, providing a means to explore 
the dynamics of novel electronic states with unprece- 
dented resolution. The possibilities include: the physics 
of the metal-insulator transition, where charge localiza- 
tion should lead to frequency scaling of the conductivity; 
electron spin resonance spectroscopy of crystal field exci- 
tations; ferromagnetic resonance in novel magnetic struc- 
tures; and cyclotron resonance studies of Fermi surface 
topology. In addition to the cuprate superconductors, 
other natural possibilities in the area of superconduc- 
tivity include heavy fermion and ruthenate materials, as 
well as high resolution spectroscopy of low-frequency col- 
lective excitations such as Josephson plasmons and order- 
parameter collective modes. 

Early microwave measurements on high quality single 
crystals of YBa2Cu307_5 showed that cooling through 
Tc « 90 K decreased the surface resistance very rapidly, 
by four orders of magnitude at 2.95 GHz, reaching a low 
temperature value of several fxft.^ Resolving this low ab- 
sorption in the microwave region has provided a technical 
challenge that has been successfully met over most of the 
temperature range below Tc by the use of high precision 
cavity-perturbation techniques'*^ In these experiments, 
the sample under test is brought into the microwave fields 
of a high quality-factor resonant structure made from 
superconducting cavities or low-loss dielectric pucks. A 
limitation of such techniques is that the resonator is gen- 
erally restricted to operation at a single fixed frequency, 
therefore requiring the use of many separate experiments 
in order to reveal a spectrum. Furthermore, a very gen- 



eral limitation of the cavity perturbation method is that 
the dissipation of the unknown sample must exceed the 
dissipation of the cavity itself in order to be measured 
with high precision — a very strong demand for a high 
quality superconductor in the T — > limit. The mea- 
surement of the residual absorption in superconductors 
is challenging at any frequency: in the case of infra-red 
spectroscopy the problem becomes that of measuring val- 
ues of reflectance that are very close to unity. The chal- 
lenge lies in the calibration of the measurement, and in 
both microwave and infra-red work, one relies on having 
a reference sample of known absorption to calibrate the 
loss in the walls of the microwave resonator or the infra- 
red reflectairce. Despite these limitations, resonant mi- 
crowave techniques are the only methods with sufficient 
sensitivity to measure the evolution of the microwave ab- 
sorption over a wide temperature range. In a recent effort 
by our group, five superconducting resonators were used 
to map a coarse conductivity spectrum from 1 GHz to 
75 GHz in exceptionally clean samples of YBa2Cu306.99 
from 4 K to 100 K.* This work revealed low temperature 
quasiparticle dynamics inconsistent with simple models 
of d-wave superconductivity, whose key signatures appear 
in the frequency dependence of the conduct ivitynS^S*'' The 
failure of simple theories to give a complete description 
of the temperature evolution and shape of the conduc- 
tivity spectra in the best quality samples has driven us 
to develop the technique described here, with the result 
that we can now resolve low-temperature microwave con- 
ductivity spectra in unprecedented detail. 

Bolometric detection is a natural method for measur- 
ing the surface resistance spectrum over a continuous fre- 
quency range. For any conductor, the power absorption 
in a microwave magnetic field is directly proportional to 
the surface resistance Rg'. 



Rs 



(1) 



where H^f is the r.m.s. magnitude of the tangential mag- 
netic field at the surface S. As a result, a measurement of 
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the temperature rise experienced by a weakly-thermally- 
anchored sample exposed to a known microwave mag- 
netic field Hrf directly gives Rs- To enhance rejection 
of spurious temperature variations, the rf power should 
be amplitude modulated at low frequency and the re- 
sulting temperature oscillations of the sample detected 
synchronously. 

We note that as part of a pioneering study of super- 
conducting Al, a similar bolometric microwave technique 
was used by Biondi and Garfunkel to examine the de- 
tailed temperature dependence of the superconducting 
gap frequencyiS This earlier experiment had the simpli- 
fying advantage of measuring the absorption by a large 
waveguide made entirely from single crystalline Al. Un- 
fortunately, in more complicated materials such as the 
multi-elemental cuprate superconductors, the best qual- 
ity samples can only be produced as small single crys- 
tals. More recently, frequency-scanned bolometric mea- 
surements have proven useful in probing collective excita- 
tions in small samples of high-Tc cuprates at frequencies 
above 20 GHz where the absorption is larger and much 
easier to measure.^ These techniques, however, have not 
focussed on the challenge of resolving the low temper- 
ature absorption of high-quality single crystals across a 
broad frequency range. 

A characteristic feature of many electronic materials 
of current interest is reduced dimensionality, which gives 
rise to highly anisotropic transport coefficients. When 
making microwave measurements, a well-defined geom- 
etry must be chosen in order to separate the individ- 
ual components of the conductivity tensor, and also to 
ensure that demagnetization effects do not obscure the 
measurement. One particularly clean approach that has 
been widely used places the sample to be characterized 
near a position of high symmetry in a microwave enclo- 
sure, in the quasi-homogeneous microwave magnetic field 
near an electric node. Often, single crystal samples grow 
naturally as platelets having a broad d~b plane crystal 
face and thin c-axis dimension, and demagnetization ef- 
fects are minimized if the broad face of the sample is 
aligned parallel to the field. In response to the applied rf 
magnetic field, screening currents flow near the surface 
of the sample along the broad a or 6 face and must nec- 
essarily flow along the c direction to complete a closed 
path. In some cases, it is desirable to work with samples 
that are very thin, rendering the c-axis contribution neg- 
ligible. Alternatively, by varying the aspect ratio of the 
sample by either cleaving or polishing, one can make a se- 
ries of measurements to disentangle the different crystal- 
lographic contributions, without having to change sam- 
ples. For example, in the cuprate superconductors, the 
conductivity parallel to the two dimensional Cu02 plane 
layers can be several orders of magnitude larger than that 
perpendicular to the weakly-coupled planes. Typical as- 
grown crystal dimensions are 1.0 x 1.0 x .01 mm^. For 
this aspect ratio, experiments where a sample was cleaved 
into many pieces showed that the c-axis contribution 
is unimportant.^- We note here that all measurements 



presented in this article employ the low-demagnetization 
sample orientation discussed above. 

The broadband surface-resistance measurement tech- 
nique we describe in the following sections provides three 
distinct technical advances over previous bolometric ap- 
proaches: a uniform microwave fleld configuration in the 
sample region that permits the separation of anisotropic 
conductivity components; the use of an in-situ reference 
sample that calibrates the microwave field strength at 
the sample absolutely; and very high sensitivity afforded 
by the choice of a resistive bolometer optimized for the 
low-temperature range and mounted on a miniaturized 
thermal stage. These features of our apparatus permit 
precision measurements of the absolute value of i?s(w, T) 
in very-low- loss samples down to 1.2 K and over the fre- 
quency range 0.1-21 GHz. We will briefly demonstrate 
that this range captures the key frequency window for 
long-lived nodal quasiparticles in extremely clean sam- 
ples of YBa2Cu307_5, and to further demonstrate the 
performance and versatility of the apparatus, we also 
show an example of zero-field electron-spin-resonance 
spectroscopy. 
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FIG. 1: A simple thermal model consisting of a heat capac- 
ity C thermally isolated from base temperature To by a weak 
thermal link of conductance K. The resistive bolometer is 
thermally anchored to C and monitors its temperature, Tb , 
which is elevated above To by a constant current bias Ib pass- 
ing through the bolometer. The absorption of incident signal 
power Ps causes heating in C, detectable as a temperature 
rise through a change in the voltage Vb- 



II. BOLOMETRIC DETECTION 

It is instructive to calculate the minimum power de- 
tectable by a simple thermal stage, the temperature 
of which is monitored by a resistive bolometer, as de- 
picted in Fig. The bolometer has a resistance Rb 
and is in thermal equilibrium with a larger heat capac- 
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ity C representing contributions from the sample, its 
holder, and the weak thermal link. This combination 
is weakly connected, via a thermal conductance K, to 
a heat sink maintained at base temperature Tq. The 
bolometer is heated to its operating temperature Tb by 
a bias power Pb = I%Rb, where Ib is the fixed bolome- 
ter bias current. For this analysis we do not consider 
feedback effects, although they are very important in 
the special case of transition edge bolometersii^ As a 
result, we consider a configuration where Ib provides 
only modest self-heating of the bolometer, such that 
7 = {Tb — To) /To < 1. An incident signal power Ps 
raises the temperature by an amount 5Tb = Ps/K, caus- 
ing a change in the readout voltage across the bolometer 
5Vb = lB{dRB/dT)STB = Ib[ARb /dT)Psl K . We then 
define a threshold detectable signal level that is equal 
to the thermal noise u„ generated in a bandwidth Ai/ 
in the bolometer, = ^UbTbRb^v. It is then possi- 
ble to write an expression for the minimum detectable 
power p™*" in terms of the dimensionless sensitivity of 
the bolometer Sd = T/i? |di?/dT|, typically of the or- 
der of unity, the noise power P„ = ksTB^v, and the 
bolometer bias power Pb'- 

Ps"' = -^VpJ^- (2) 

From this expression one immediately sees that it is desir- 
able to minimize both the bias and noise powers, within 
the combined constraints of maintaining the bolometer 
temperature at Tb and keeping the thermal response 
time fixed at a suitably short value. By miniaturization 
of the sample holder, the bias power required to reach 
a given temperature can be considerably reduced, while 
at the same time maintaining a practical thermal time 
constant. The noise power is limited intrinsically by the 
thermal (Johnson) noise from the bolometer resistance at 
temperature Tb- However, most real sensors show sub- 
stantial excess noise, and the Cernox 1050 sensor-'^'* used 
in the present implementation of our experiment is no 
exception, showing approximately 40 dB of excess noise 
in the presence of a 1.3 /xA bias current. This completely 
accounts for the discrepancy between the minimum de- 
tectable power at 1.3 K of 17 fW calculated using Eq. [21 
assuming only Johnson noise, and the experimentally de- 
termined value of 1.5 pW. 

III. BROADBAND ABSORPTION 
MEASUREMENT APPARATUS 

For our method of bolometric detection to be most 
useful, it is necessary to deliver microwaves to the sam- 
ple across a broad range of frequency and, at the same 
time, not only accurately control the polarization of the 
microwave field at the sample, but also maintain a fixed 
relationship between the field intensity at the sample un- 
der test and the field intensity at the reference sample. 
Essential to this is the design of the microwave waveg- 



uide. We use a custom-made transmission line, shown 
in cross-section in Fig. |21 that consists of a rectangu- 
lar outer conductor that measures 8.90 mm x 4.06 mm 
in cross-section and a broad, flat centre conductor, or 
septum, that measures 4.95 mm x 0.91 mm. This sup- 
ports a TEM mode in which the magnetic fields lie in the 
transverse plane and form closed loops around the centre 
conductor, setting a fixed relationship between the mi- 
crowave field strengths on either side of the septum. The 
line is terminated by shorting the centre conductor and 
outer conductor with a flat, metallic endwall. This en- 
forces an electric field node at the end of the waveguide, 
adjacent to which we locate the small platelet sample and 
reference, with their flat faces parallel and very close to 
the endwall. The broad centre conductor ensures spa- 
tially uniform fields over the dimensions of the sample, 
making it possible to drive screening currents selectively 
along a chosen crystallographic direction. The electro- 
dynamics of the rectangular waveguide are discussed in 
more detail in Appendix ^ A strong variation in the 
power delivered to the sample as a function of frequency 
arises due to standing waves in the microwave circuit. 
In order to properly account for this, we have incorpo- 
rated an in-situ normal-metal reference sample of known 
surface resistance that acts as an absolute power meter. 
This second sample is held in a position that is electro- 
magnetically equivalent to the that of the test sample, 
on a separate thermal stage. 

copper Cernox sapphire sample 

holder sensor plate under test 

d . L _i_ I 

t t I 7 ■ 

quartz chip reference 

tube heater jgrnple 

FIG. 2: Schematic cross-section of the terminated coaxial line 
region showing the sample and reference materials suspended 
on sapphire plates in symmetric locations in the rf magnetic 
field. The sapphire plate is epoxied into the bore of a quartz 
tube which thermally isolates it from the copper holder, fixed 
at the temperature of the 1.2 K helium bath. 

One of the challenges of cryogenic microwave absorp- 
tion measurements on small, low-loss samples is the de- 
sign of the sample holder, which must measure and reg- 
ulate the sample temperature, and yet contribute neg- 
ligible dissipation compared to the sample. A widely 
used technique that satisfies these requirements is that 
of a sapphire hot-finger in vacuum, allowing the ther- 
mometry to be electromagnetically shielded from the mi- 
crowave fields. In our apparatus, the sample holder is 
inserted through a hole that is beyond cut off for all op- 
erating frequencies. For ac calorimetric measurements, 
the design of the thermal stage is critical in setting the 
sensitivity of the system. The experimental arrangement 
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RF coaxial line 

copper vacuum can 
stainless steel tube 
Cernox sensor 

nylon spacer 
sapphire plate 

HITc sample 
center conductor 
sapphire plate 

chip heater 
Cernox sensor 
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indium seal 



FIG. 3: Scale drawing of the assembled apparatus indicating 
the details of the vacuum can and sample region. The alloy 
reference sample is not visible in this cut-away view. 



is shown schematically in Fig. [21 with the sample under 
test fixed on the end of a 100 /im thick sapphire plate us- 
ing a tiny amount of vacuum greased The plate extends 
17 mm from the sample to where it is epoxied into the 
bore of a 1.2 mm diameter quartz glass tube that acts as 
a thermal weak-link to the liquid helium bath. A Cernox 
thermometer and a 1500 surface-mount resistor used 
as a heater— are glued directly onto the sapphire plate 
with a very thin layer of Stycast 1266 epoxy,^^ ensuring 
intimate thermal contact with the sapphire and hence the 
sample. We use 0.05 mm diameter NbTi superconduct- 
ing electrical leads to the thermometer and heater for 
their very low thermal conductance, which is in parallel 
with the quartz weak-link. 

The microwave circuit is powered by a Hewlett- 
Packard 83630A synthesized sweeper (0.01-26.5 GHz) 
combined with either an 8347A (0.01-3 GHz) or 8349B (2- 
20 GHz) amplifier, generating up to 23 dBm of rf power 
across the spectrum. Approximately 2 m of 0.141" stain- 
less steel coaxial linei^ delivers power from the amplifier 
down the cryostat to the vacuum can where it is soldered 
into the rectangular line (see Fig. The r.m.s. mi- 
crowave magnetic field amplitude at the samples is typ- 
ically ~10~^ Oersteds, which generates ~ /iK modula- 
tions in the sample-stage temperature for a typical high 
quality 1 mm^ high-Tc sample having a low frequency Rs 
value of 1 fifl. 

An assembled view of the low temperature apparatus 
including the microwave transmission line and the posi- 
tions of the sample and reference holders is provided in 



Fig. O The sapphire plates that support both the test 
sample and reference sample are inserted through 4 mm 
cut-off holes into the microwave magnetic field. The rect- 
angular coaxial line consists of a centre conductor made 
from a 0.91 mm thick copper plate soldered at one end 
onto the centre conductor of the 0.141" semi-rigid coax- 
ial line, and at the other end into the wall of the copper 
cavity that comprises the outer conductor of the trans- 
mission line. To minimize the rf power dissipated in the 
low temperature section of transmission line, the entire 
surface exposed to microwave radiation, including the fi- 
nal 15 cm of semi-rigid coaxial line, was coated with PbSn 
solder, which is superconducting below 7 K. During ex- 
periments, the vacuum can is completely immersed in a 
pumped liquid helium bath having a base temperature of 
1.2 K. 

The selection of a reference material for low-frequency 
work must be made carefully. Initially, we chose sam- 
ples cut from commercially available stainless-steel shim 
stock, a common choice in infrared spectroscopy work. 
Calibration experiments produced erratic results which 
were eventually traced to the presence of anisotropic 
residual magnetism in the stainless steel. Subsequently, 
we produced our own reference material, choosing an 
Ag:Au alloy (70:30 at.% made from 99.99% pure starting 
materials), because it exhibits a very simple phase dia- 
gram that guarantees homogeneity— By using an alloy, 
we ensure that the electrodynamics remain local at mi- 
crowave frequencies, avoiding the potential complications 
arising from the anomalous skin effect.^*' Our sample was 
cut from a 93±5 /im thick foil having a measured residual 
dc resistivity value of /9=5.28±0.3 /ii7cm, constant below 
20 K. 

While the thermal stage for the reference sample is 
similar in design to that used for the sample under test, 
it uses a higher conductance stainless steel thermal weak- 
link (in place of the quartz tube), since the dissipation of 
the normal metal calibration sample is orders of magni- 
tude larger than that of a typical superconducting sam- 
ple. Because the apparatus was implemented as a retro- 
fit to an existing experiment, the reference thermal stage 
had to be mounted directly onto the body of the trans- 
mission line structure. Although the cavity walls are su- 
perconducting to reduce their absorption, we use a nylon 
spacer to thermally isolate the base of the reference from 
the transmission line to avoid direct heating. The heat- 
sinking of the reference base to the helium bath is made 
using a separate copper braid that is not visible in Fig.|31 

As considered previously in our generic analysis, we op- 
erate the Cernox bolometer with a constant dc current 
bias, typically a few fiA, provided by the series combi- 
nation of an alkaline battery (1.5 V or 9 V) and bias 
resistor whose value is much larger than that of the Cer- 
nox sensor. With the helium bath under temperature 
regulation, the choice of bias power sets the tempera- 
ture of the sample for a given experiment, with no other 
temperature control necessary. All electrical leads into 
the cryostat are shielded twisted pairs of insulated man- 
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ganin wire, and true four-point resistance measurements 
are made on all sensors. The voltage signal appearing on 
the Cernox thermometer is amplified outside the cryo- 
stat by a carefully shielded and battery-powered circuit. 
We use a two-stage cascaded amplifier with one Analog 
Devices AD548 operational amplifier per stage, chosen 
because these are readily available, low-noise amplifiers. 
The dc level is nulled between stages to prevent satura- 
tion, and the total gain is 10*. The amplified signal, cor- 
responding to the temperature modulation of the sample, 
is then demodulated with a Stanford Research Systems 
SR850 digital lock-in amplifier that is phase-locked with 
the rf-power amplitude modulation. There are two such 
systems, one for the sample and one for the reference 
measurements. The entire experiment is operated under 
computer control when collecting data. 



IV. CALIBRATION 

Two steps are necessary for an absolute calibration 
of the surface resistance of an unknown specimen from 
the measured temperature-rise data. The first is to cal- 
ibrate the absolute power sensitivity of the sample and 
reference thermal stages at the actual operating temper- 
ature and modulation frequency. This is achieved using 
the small in-situ heater to drive well-characterized heat 
pulses that mimic absorption by the sample, while at the 
same time measuring the corresponding temperature re- 
sponse. The second step requires the calibration of the 
magnetic field strength at the sample, at each frequency, 
using the known absorption of the reference sample. We 
exploit the fact that the metallic reference sample experi- 
ences the same incident microwave field Hrf as the sam- 
ple under test, guaranteed by conservation of magnetic 
flux. This ensures that taking the ratio of the absorbed 
power per unit surface area of each sample provides the 
ratio of the surface resistance values: 

rjsam rysam A sam 

pref U^'^f Aref ' ^ ' 

The surface resistance of the unknown sample i?™'"(a;) is 
then trivially found by multiplying the power-absorption 
ratio, shown in Fig. ^ by Rg (w) of the metallic reference 
sample cal culated using the classical skin-effect formula 
Rs{uj) = ujfiop/2 where Lojlix is the frequency and /io is 
the permeability of free space. The raw power-absorption 
spectra, shown in the first two panels of Fig. 0] highlight 
the necessity of the reference sample. The absorption 
spectra of the samples is completely masked by the large 
amplitude variations of H^f caused by the standing waves 
in the microwave circuit. 

An essential test of the method is to make a frequency- 
scanned measurement with identical samples mounted on 
the sample and the reference stages. The result should 
be a frequency-independent ratio across the spectrum, 
equal to unity for samples with the same surface area. 
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FIG. 4: Raw absorption spectra corresponding to the tem- 
perature rise of the sample. Taking the ratio of the two sig- 
nals accounts for the strong frequency dependence of Hrj in- 
troduced by standing waves in the transmission line. The 
remaining frequency dependence of the ratio is due to the 
different Rs(y^) spectra of the two samples. 



In Fig. [S] we show such a measurement using two thin- 
platelet samples of our Ag:Au reference alloy, both at 
the base temperature of 1.2 K. The data reveal a ratio of 
0.82 in the cryogenic apparatus, due to the fact that the 
centre conductor is offset from centre by 0.1 mm in the 
termination region, intensifying the fields on one side rel- 
ative to the other. This scale factor must be included in 
the calibration of all experimental data. The sharp reso- 
nance seen in the ratio at 22.5 GHz indicates the presence 
of a non-TEM electromagnetic mode in the sample cav- 
ity that breaks the symmetry in field strength between 
sample and reference positions. For our present design. 
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FIG. 5: Ratio of the sample absorption to reference absorp- 
tion for identical samples, compared to measurements of the 
field amplitude at equivalent positions in a 4x scale model. 
(Frequencies for the scale model have been scaled by a fac- 
tor of 4 for the comparison.) The ratio technique is seen to 
break down with a sharp resonance in both cases. The origin 
of these resonances, which limit the useful frequency range 
of the apparatus, is discussed in detail in Appendix ^ and 
shown to be due to the presence of standing waves of the TEoi 
waveguide mode. For this mode, the three arrows indicate: 
the cut-off frequency f^^ — 15.38 GHz at which the mode is 
first free to propagate; its quarter-wave resonance frequency 
/^y4 = 19.11 GHz, for open-circuit termination conditions; 
and the half-wave resonance frequency /^'jij — 27.4 GHz, for 
short-circuit termination conditions. The A/4 and A/2 reso- 
nance frequencies bracket the observed resonances. The scale 
model, which has a large transition capacitance between cir- 
cular and rectangular coax sections, is seen to fall at the high 
end of the range. 



this sets the upper frequency limit of operation. 

In an attempt to gain further insight into the field con- 
figurations in the transmission hne, and to understand 
how the higher order waveguide modes hmit the upper 
frequency range, we built a scale model of the setup hav- 
ing all dimensions larger than those of the cryogenic ap- 
paratus by a factor of four. For comparison, a frequency 
scan of the model structure is included in Fig. [S] us- 
ing loop-probes in the positions of the samples. The 
data show that the non-TEM-mode resonance occurs at 
27 GHz, considerably higher than in the low tempera- 
ture experiment. It turns out that the breakdown of the 
sample-reference symmetry occurs not at the frequency 
at which higher order waveguide modes first propagate in 
our structure, but at the frequency at which they form a 
resonant standing wave. A full discussion of this is given 
in Appendix IbI 

A number of other experimental tests were important 
to verify the proper operation of the system. Frequency 
scans without samples mounted on the sapphire stages 
confirmed that background absorption due to the sap- 
phire and tiny amount of vacuum grease used to affix 
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FIG. 6: Comparison of measurements made on the same 
sample of YBa2Cu3O6.50 using the broadband bolometric ex- 
periment (solid symbols) with those from five microwave res- 
onators (open symbols). The agreement between methods is 
excellent. The data is plotted as Rb{uj)/lo'^ to remove the 
frequency dependence associated with superfluid screening. 



the samples is negligible — it is unmeasurable at low 
frequency, and contributes no more than 2 fifl to an Rs 
measurement at 21 GHz. Scans without a sample also 
confirmed that no significant leakage heat current prop- 
agates to the thermometers directly from the microwave 
waveguide. The high thermal stability of the cryostat 
system is due in part to the very large effective heat ca- 
pacity of the pumped 4 litre liquid-helium bath at 1.2 K. 
In addition, it is always important to make certain the 
temperature modulations of the samples are sufficiently 
small that the response of the thermal stages remains in 
the linear regime. Furthermore, measurements with the 
same sample located in different positions along the sap- 
phire plate, with up to 0.5 mm displacement from the 
central location in the waveguide, confirmed that there 
is enough field homogeneity that our sample alignment 
procedure using an optical microscope is sufficient, and 
that samples of different sizes experience the same fields. 

A very convincing verification of the technique is pro- 
vided by the ability to compare broadband Rs{iO-, T) data 
with measurements of the same sample in five differ- 
ent high-Q microwave resonators. These experiments 
probe the temperature dependence of the absorption to 
high precision at a fixed microwave frequency: however, 
the determination of the absolute value of Rs is limited 
to about 10% as discussed previously. The bolometric 
method has the advantage of being able to measure a 
true spectrum because the dominant uncertainty, the ab- 
solute surface resistance of the reference sample, enters 
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as a scale factor that applies across the entire frequency 
and temperature range. A detailed discussion of the un- 
certainties in the bolometry data will be presented in the 
subsequent section. Figure El shows that there is very 
good agreement of both the temperature and frequency 
dependence of the surface resistance as measured inde- 
pendently by the fixed-frequency and broadband experi- 
ments. 
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FIG. 7: Broadband measurements of the microwave surface 
resistance spectrum of YBa2Cu3O6.50 obtained with the bolo- 
metric apparatus below 10 K. The low frequency absorption 
approaches the resolution limit of the apparatus, while the 
upper frequency limit is imposed by the resonance in the mi- 
crowave structure. 



V. PERFORMANCE 

Figure [71 presents an example of high resolution broad- 
band measurements of the frequency-dependent and 
temperature-dependent surface resistance of a supercon- 
ducting sample. This particular data set is for d-axis 
currents in a YBa2Cu306.5o single crystal (Tc=56 K) 
having dimensions 1.25x0.96x0.010 mm"^. The data 
span the range 0.6-21 GHz, hmited at high frequency 
by the resonance in the system, and at low frequency by 
the small dissipation of the sample, which approaches the 
resolution limit of the experiment. At 1 GHz, the values 
for the statistical r.m.s. uncertainty in surface resistance, 
6Rs, are about 0.2, 0.4, 0.6, and 1.3 ^mQ for T = 1.3, 2.7, 
4.3, and 6.7 K respectively. Error bars have been omit- 
ted from the figure for clarity. Systematic contributions 
to the uncertainty enter as overall scale factors in the Rs 
data and are attributed to an uncertainty in the DC re- 
sistivity of the thin Ag:Au alloy foil used as a reference 
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Frequency (GHz) 

FIG. 8: The real part of the microwave conductivity (Ji{uj, T) 
extracted from the broadband Rs{to,T) measurements. We 
use the self-consistent fitting procedure described in Ap- 
pendix |ni to properly account for contributions to field- 
screening by the quasiparticles. 



sample (~5%), the surface area of the samples (~1%), 
and the absolute power sensitivity of the thermal stage 
(-1%). 

The frequency dependence observed in Rs{w) is due 
to absorption by quasiparticles thermally excited from 
the superfiuid condensate. The quantity of fundamental 
theoretical interest is the real part of the conductivity 
spectrum ai (cj), which must be extracted from the exper- 
imentally measured Rs{uj) data. A thorough discussion 
of the method we use to do this is given in Appendix lO 
but, to first approximation, the shape of the conductivity 
spectrum can be found by dividing by a factor of 

to account for the screening of the applied field by 
the superfiuid. Figure |H1 shows the conductivity spectra 
extracted from the Rs{uj) data using the complete anal- 
ysis. It is immediately apparent why improving the sen- 
sitivity of the experiment is of the utmost importance. 
The low frequency region, where the power absorption 
becomes very small, is where the conductivity exhibits 
the strongest frequency dependence and is most impor- 
tant to measure accurately. The spectrum at 1.3 K has a 
width of the order of 5 GHz, signifying very long quasi- 
particle scattering times, indicative of the high quality of 
our YBa2Cu306.5o crystal. For this very clean sample, 
most of the spectral weight resides below the experimen- 
tal frequency limit of 21 GHz at 1.3 K, but the increase 
in scattering with increasing temperature quickly broad- 
ens the spectra thus motivating future designs capable of 
probing a broader frequency range. Many cuprate ma- 
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FIG. 9: The electromagnetic absorption spectrum of a 
Y0.99Gd0.01Ba2Cu3O6.99 single crystal at 1.3 K. The spec- 
trum consists of a broad background due to the quasiparticle 
absorption in the superconductor in addition to the zero field 
ESR lines generated by the low concentration of magnetic 
Gd'^^ impurities. Only ions residing within a distance A of 
the crystal surface contribute to the signal as the applied rf 
field is strongly screened by the superconductor. 



points that occur in certain d- and /-electron metals. 
With a httle attention to thermal design, specifically the 
thermal separation of sample stages from the microwave 
waveguide, the superconducting coatings on the waveg- 
uide could be removed and the system used in high mag- 
netic fields. This would open interesting possibilities in 
the area of metals physics, such as high-resolution cy- 
clotron and periodic-orbit resonance, as well as the study 
of vortex dynamics and vortex-core spectroscopy in su- 
perconductors. Finally, further miniaturization of the 
experiment should also be possible: the ultimate goal 
would be to extend the frequency range of this type of 
spectroscopy so that is joins seamlessly on to the THz 
range now accessible using pulsed-laser techniques. 
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terials, such as Bi2Sr2CaCu208+5, have scattering rates 
that are orders of magnitude higher and require THz fre- 
quency techniques to probe their dynamics. 

As a final demonstration of the sensitivity of the broad- 
band instrument we have described, we include a fre- 
quency scan of a superconducting sample that exhibits 
clearly discernable absorption lines due to zero field elec- 
tron spin resonance (ESR) of a low density of magnetic 
impurities. Gadolinium ions (electron spin S = 7/2) sub- 
stitute for yttrium, sandwiched between the two CuO 
planes in the YBaCuO unit cell, and the splitting of 
the degenerate Gd^+ hyperfine levels by the crystalline 
field provides a very sensitive probe of the local micro- 
scopic structured These measurements are typically per- 
formed in a high field spectrometer, but the bolometric 
system provides a means of performing zero-field mea- 
surements. Figure El shows the 1.3 K absorption spec- 
trum of a 1 mm^ Y0.99Gd0.01Ba2Cu3O6.99 single crystal 
consisting of a broad background due to the quasipar- 
ticle conductivity, essentially unaltered by the presence 
of the Gd ions, with the ESR absorption peaks super- 
posed. The high signal-to-noise ratio achieved with the 
experiment allows one to resolve the spectrum in great 
detail. 

The apparatus described here has sufhcient sensitivity 
and frequency range for it to be immediately applicable 
to many other interesting problems in condensed matter 
physics. These might include: the study of low- lying col- 
lective modes in metals and superconductors; zero-field 
electron spin resonance in correlated insulators; and the 
study of critical phenomena at the metal-insulator tran- 
sition and near the zero-temperature magnetic critical 



APPENDIX A: FREQUENCY RESPONSE OF 
DISTRIBUTED THERMAL STAGE 

We wish to calculate the temperature response of a 
simple thermal stage to a sinusoidal heat flux Pi„ — 

RejPoe^"*} superimposed on a static temperature gra- 
dient. We consider an arrangement where the isother- 
mal sample stage has neglible thermal mass and is con- 
nected to base temperature by a weak thermal link 
with distributed heat capacity cy per unit volume. It 
is a straightforward extension to include an additional 
lumped heat capacity for the isothermal stage; once the 
lumped heat capacity dominates, the frequency response 
simplifies to that of a single-pole low-pass filter. How- 
ever, in our case this is unnecessary: for electrodynamic 
measurements at low temperatures, the sample holder 
is required to be both electrically insulating and highly 
crystalline, and will therefore have very low heat capac- 
ity. 

Here we consider the one-dimensional problem of a 
thin bar (the quartz tube in our apparatus) of length 
£ and cross-sectional area A, with one end fixed at a base 
temperature To while the other end is heated by a heat 
flux due to sample power absorption. The propagation 
of a heat current Jq through the bar is constrained by 
the continuity equation d-Jq/dx + cydT/dt = and the 
thermal conductivity k is defined by Jq = —ndT / dx. 
Together, these lead to the one-dimensional heat equa- 
tion dT/dt = ad^T/dx^ where a = k/cv is the thermal 
diffusivity. Defining a complex thermal diffusion length 
5 = ^J a/iw, the time-dependent part of the temperature 
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profile can be written 

AT{x, t) = Re jfe^"* sinh(x/^)| 



(Al) 



where T is fixed by the heat-flux boundary condition: 
Pin = —KA{dT/dx)\x=i- This completely determines 
the frequency-dependent temperature rise of the sample 
stage: 



AT(i, i) = Re i ^Se''^' tanh(£/(5) 
kA 



1.1 



(A2) 



In the low frequency limit, the temperature rise reverts 
to the usual result: 



Ar(£,t) = ^cosM, 



(A3) 



where, without loss of generality, we have set the phase 
of the input heat flux to zero. In the high frequency 
limit, the thermal diffusion length becomes shorter than 
the weak link and the temperature rise is reduced, being 
given by: 



ATi£, t) = cos{ujt - 7r/4) 



(A4) 



where \6\ = y^a/uj. At finite frequencies, part of the 
heat flux is diverted into the distributed heat capacity 
of the thermal link. For a fixed input power (and hence 
fixed temperature gradient at the end of the thermal link) 
this leads to smaller temperature rises and a decreased 
sensitivity of the bolometric method. Clearly, the ex- 
perimental sensitivity of the bolometric method will be 
optimized by operating in the low- frequency limit: > £ 
or oj < a/£'^. A consideration of the thermal diffusivity 
and dimensions of the weak-link must therefore be part 
of any plan to increase the modulation frequency. 

Fig. ^1 shows the frequency response of the sam- 
ple thermal stage in our apparatus when it was sub- 
jected to a sinusoidally varying heater power, normal- 
ized to the static response. Included in the figure are 
fits to the distributed-heat-capacity model, Eg . IA2I and a 
single-pole low-pass filter response AT{uj)/AT{uj = 0) = 
1 / Vl + uj'^t'^- Although both curves fit the data well over 
most of the frequency range, the best-fit value of the time 
constant r in the lumped-element model corresponds to 
a heat capacity much larger than the calculated heat ca- 
pacity of the sapphire sample stage. Instead, the value 
obtained from the fit is approximately half the heat ca- 
pacity of the quartz tube, indicating the correct physics 
is that of heat diffusion in a distributed thermal system. 



APPENDIX B: DESIGN STRATEGY FOR 
RECTANGULAR COAXIAL TRANSMISSION 
LINE 

In optimizing a microwave transmission line for the 
bolometric measurement of surface resistance the guid- 
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FIG. 10: Low temperature (T = 1.3 K) measurements of 
the dynamic thermal response of the quartz-tube bolometer 
platform. Curves on the plot show fits using lumped and 
distributed heat capacity models. 



ing aims must be: to deliver microwave power effi- 
ciently to the sample region, over as wide a frequency 
range as possible and with a well defined polarization; 
to have regions of uniform microwave magnetic field at 
the sample and reference positions; and, at these posi- 
tions, to have a fixed, frequency-independent ratio be- 
tween the field strengths. These aims can be met by us- 
ing an impedance-matched (50 Q), single-mode coaxial 
line, with rectangular cross section and a broad, flat cen- 
ter conductor or septum, as shown in Fig. lllf il. In addi- 
tion, the dimensions of the rectangular coaxial line should 
be chosen carefully, to prevent higher-order waveguide 
modes from entering the operating frequency range of the 
experiment, as these modes break the symmetry in field 
strength between sample and reference positions. This 
appendix outlines how to undertake the optimization. 

A rectangular coaxial line, like any two-conductor line, 
supports a transverse electromagnetic (TEM) wave at all 
frequencies. Figure lTTT ii) shows its electric and magnetic 
field configurations. The TEM mode has the desirable 
property that its magnetic fields lie in a plane perpen- 
dicular to the direction of propagation, forming closed 
loops around the centre conductor. Conservation of mag- 
netic fiux then leads to a fixed, frequency-independent 
relation between the fields on either side of the septum. 
These fields will also be quite homogeneous, as long as 
the height C of the centre conductor is large compared 
to the gap {B — D) /2 between the centre and outer con- 
ductors. To deliver microwave power efficiently to the 
sample region the characteristic impedance of the TEM 
mode must be close to that of the cylindrical coaxial line 
used to bring microwaves into the cryostat. Gunston^ 
has tabulated data on the impedance of rectangular coax- 
ial line, and gives some useful approximate formulas. The 
following expression, due to Brackelmann, is stated to be 
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TEM 




► Lines of electric flux 

► Lines of magnetic flux 

FIG. 11: Cross sections of the rectangular coaxial transmis- 
sion line: (i) physical layout of the transmission line show- 
ing dimensions {A, B ,C, D) of the inner and outer conduc- 
tors, and sample and reference positions (shaded squares) ; (ii) 
fields of the TEM mode, showing how continuity of flux and 
a broad, flat inner conductor produce uniform, well polarized 
fields of equal intensity and opposite direction at the sam- 
ple and reference positions; (iii) fields of the TEio mode; and 
(iv) fields of the TEoi mode. (Magnetic fields of the trans- 
verse electric modes contain a component along the direction 
of propagation and do not form closed loops in the transverse 
plane.) It is clear that the TEoi mode is most harmful to the 
operation of the broadband apparatus: its magnetic fields 
have high intensity at the sample and reference positions and 
break the balance that otherwise exists in the TEM mode. 



accurate to 10% for D/B < 0.3 and C/A < 0.8: 

ZoV^ = 59.952 In (^^i^^f], (Bl) 

where is the relative permittivity of the dielectric filling 
the transmission line. 

We now come to the question of what places an up- 
per limit on the useful frequency range of the rectangu- 
lar coaxial waveguide. At high frequencies our method, 
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FIG. 12: Conducting walls introduced along special electric 
equipotentials allow the waveguide modes of rectangular coax- 
ial line to be mapped onto the fundamental mode of ridged 
waveguide, a problem extensively studied in the literature. 
The figures show the relabelling of dimensions in Pyle's nota- 
tion-S, as a, 6, d and s. 



which incorporates an in-situ power meter, suffers a spec- 
tacular breakdown in the ratio of the relative strengths 
of the microwave magnetic fields at the sample and ref- 
erence positions, as shown in Figure |S1 This is caused 
by the presence of higher-order waveguide modes, which 
have different character from that of the TEM mode un- 
der the mid-plane reflection symmetries of the rectangu- 
lar line. The waveguide modes with the lowest cut-off 
frequencies are the transverse electric modes TEio Sind 
TEoi , shown in Figures lllf iii) and lllf iv) respectively. 
These have the property that magnetic fields on oppo- 
site sides of the septum point in the same direction. The 
fields of the TEM mode, in contrast, are antiparallel, 
causing an admixture of TEM and TE modes to lack 
the important characteristic of equal field intensities at 
sample and reference positions. Particularly damaging is 
the TEoi mode, which is not screened by the septum and 
has high field intensity in the vicinity of the sample and 
reference. In principle it is possible to avoid exciting the 
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transverse electric modes by building a very symmetric 
transmission line. In practice, however, we find this to be 
impossible — sufficiently large symmetry-breaking per- 
turbations are always present. Nevertheless, maintaining 
high symmetry is still desirable. A comparison of our 
results with calculations of the cut-off frequencies of the 
transverse electric modes shows that at frequencies where 
the higher order modes are free to propagate, they do not 
immediately cause a breakdown in field ratio: this only 
occurs when the transverse electric modes come into res- 
onance. (This can be seen very clearly in Figure|Sl) As a 
result, the range of operating frequency can be extended 
by as much as 50% just by shortening the final section of 
transmission line and carefully designing the transition 
between the cylindrical and rectangular sections. 

Optimizing the range of single-mode operation of the 
rectangular transmission line requires a method for cal- 
culating the cut-off frequencies of the TEio and TEqi 
modes. While waveguide modes in two-conductor rectan- 
gular transmission lines have not been extensively stud- 
ied, their field configurations can be mapped onto a more 
common geometry: that of ridged waveguide. Figure IT^ 
shows how. Electric equipotentials run perpendicular to 
lines of electric flux, and special equipotentials, corre- 
sponding to local minima of the magnetic flux density, 
exist on the symmetry axes of the rectangle. A con- 
ducting wall can be introduced along these lines without 
disturbing the field distributions, thereby mapping each 
mode onto an equivalent ridged waveguide. Figure^] il- 
lustrates the two different ways this is done, for the TEio 
and TEoi modes respectively. A very early calculation 
of the cutoff frequency of ridged waveguide was carried 
out by Pyle2^ and is notable for its simplicity, general- 
ity and enduring accuracy when compared to more re- 
cent numerical methodsi^ Pyle's approach is to solve 
for the transverse resonance condition of the waveguide, 
which is equivalent to finding the cut-off frequency Uc- 
We have used this method in our design process, as it is 
easy to implement (involving only algebraic equations) 
and is accurate to several percent except when the sep- 
tum becomes very thin. The length £ of the rectangular 
line, the cut-off frequency, and the discontinuity capaci- 
tance of the cylindrical-to-rectangular transition together 
determine the resonant frequencies lor of the transverse 
electric modes. There are two limiting cases, correspond- 
ing to open-circuit {£ = A/4) and short-circuit {i = A/2) 
termination (where A = 27r//c is the wavelength along 
the guide) , that follow from the waveguide dispersion re- 
lation. 

2 2,2 7 2 2 I " ^oo\ 

= ujI + : open circuit (B3) 

^ ujI + : short circuit (B4) 

A high capacitance for the TE modes at the transition 
from cylindrical to rectangular coax is clearly favourable: 



it better approximates the short circuit termination con- 
dition and leads to resonant frequencies at the upper end 
of the range. This effect is responsible for the difference 
in resonant frequencies between the scale model and the 
actual apparatus seen in Figure [5| There is, however, 
a trade-off to be made: too large a transition capaci- 
tance for the TEM mode will result in most of the mi- 
crowave power being reflected before it reaches the sam- 
ple. The dimensions of the rectangular guide in our ap- 
paratus are A — 8.90 mm, B = 4.06 mm, C = 4.95 mm, 
D ~ 0.91 mm and I = 6.60 mm. The cut-off frequen- 
cies for the TEio and TEqi modes are calculated to 
be 19.68 GHz and 15.38 GHz respectively. The quar- 
ter wave- resonances would then occur at 22.72 GHz and 
19.12 GHz, and the half wave resonances at 30.06 GHz 
and 27.44 GHz. 



APPENDIX C: EXTRACTION OF a^iiS) FROM 
ii,(t^) MEASUREMENTS 

In this appendix we show how the microwave con- 
ductivity spectrum of a superconductor can be obtained 
from a measurement of its frequency-dependent surface 
resistance. This process is similar to the extraction of 
conductivity spectra in the infra-red frequency range 
from reflectance measurements. In both cases, we be- 
gin with incomplete information about the electrody- 
namic response: the bolometric technique described in 
this paper measures only the resistive part of the surface 
impedance; and optical techniques typically obtain the 
magnitude, but not the phase, of the reflectance. How- 
ever, the conductivity a = ai — ia2 is a causal response 
function, and its real and imaginary parts are related by 
a Kramers-Kronig transform: 

a.iu^) = -v r^^dn, (CI) 

where V denotes the principle part of the integral. At 
first sight we seem to have replaced one uncertainty, in- 
complete knowledge of the phase, by another, the finite 
frequency range over which the measurements have been 
made. However, a suitable extrapolation of the data out 
of the measured frequency range is usually possible and 
makes the transform a well-defined procedure in practice. 

We consider the limit of local electrodynamics, in 
which the microwave surface impedance Zg is related to 
the complex conductivity in a straightforward manner by 
the expression 

= Rs + iX, = IJ^. (C2) 

V (Tl - 1(T2 

Very generally, the conductivity can be partitioned into 
a superfluid part as, consisting of a zero- frequency delta 
function and an associated reactive term, and a normal- 
fluid component a^: 

a{u, T) = (Jis - 1(725 + 0-17V - io-2Ar (C3) 



12 



«se' .Use- 

IT d(U!) — 1 



- CTiAT — ia2N, 



where m* is the quasiparticle effective mass. In the clean- 
Umit, where the quasiparticle scattering rate l/r is much 
less than the spectroscopic gap 2A, sum-rule arguments 
enable a clean partitioning of the conduction electron 
density n into a superfluid density rig and a normal-fluid 
density Un = n — Us- In this type of generalized two-fluid 
modelfSI the temperature dependence of Ug is determined 
phenomenologically from measurements of the magnetic 
penetration depth A through the relation 



(C4) 



Applying this, we can write the conductivity at finite 
frequencies as 



(t{uo,T) = criAr(tJ,T) 



a2N{'^-,T) ■ 



1 



(C5) 

From Ea. lC2l it is clear that Rs{i^) is determined by both 
the real and imaginary parts of the conductivity. How- 
ever, one simplification occurs at temperatures well be- 
low Tc, where few thermally excited quasiparticles exist, 
and the low frequency reactive response is dominated by 
the superfluid. In this case, a good approximation to the 
relations becomes 



1 



Xs{oj,T) = hqujX{T). 



(C6) 



At higher frequencies and temperatures, a more complete 
treatment would account for quasiparticle contributions 
to field screening, which enter through tT2Jv(w,T). We 
use an iterative procedure to obtain the quasiparticle con- 
ductivity spectrum (Tijv(^), starting from the good initial 
guess provided by Eq. IC6I The process goes as follows. 
A phenomenological form that captures the key charac- 
teristics of the dataset but has no physical motivation, 
namely cri(w) = cro/[l + ("^/r)^], is fitted to the spec- 
trum and used to extrapolate out of the measured fre- 
quency range. The Kramers-Kronig transform (Eq. IC1|) 
can then be applied to obtain (T2Jv(w). With (T2Jv('jj) in 



hand, and with the superfluid contribution to a2 known 
from measurements of the magnetic penetration depth, a 
new extraction of the conductivity from the Rs{oj) data is 
made, this time using the exact expression, Eq. IC2I The 
whole procedure is repeated to self-consistency. We find 
that the procedure is stable and converges rapidly, and 
is not sensitive to the details of the high-frequency ex- 
trapolation. Also, the corrections are quite small for the 
low temperature dataset shown in Fig. [S] at the highest 
temperature and frequency they amount to a 7% change 
in ai. 

In addition, two independent experimental checks give 
us further assurance that we obtain the correct conduc- 
tivity spectra. We first note that Eq. IC2I contains an 
expression for the surface reactance Xg. Therefore, once 
we have measured the penetration depth and extracted 
the conductivity spectra from the Rs{w) data, we can 
predict the temperature dependence of the surface re- 
actance at any frequency and compare with experiment. 
We have made this comparison at 22.7 GHz, with surface 
reactance data obtained on the same crystal, and find ex- 
cellent agreement. We note that this is a frequency high 
enough for quasiparticle scattering to have a discernible 
effect on the surface reactance. 

A second verification of the conductivity extraction 
procedure is its ability to predict the spectral weight that 
resides outside the frequency window of the measure- 
ment. A corollary of the Kramers-Kronig relation ICll 
is the oscillator-strength sum rule 
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<7i{uj, T)du!. 



(C7) 



In a superconductor, in the clean limit, the sum rule re- 
quires that any spectral weight disappearing from the 
superfluid density ns(T) as temperature is raised must 
reappear as an increase in the frequency-integrated quasi- 
particle conductivity. We have carried out this compar- 
ison, which is shown in the inset of Fig. |HI^^ The good 
agreement in the temperature dependence of the super- 
fluid and normal-fluid densities is a strong veriflcation 
of both the analysis procedure and the bolometric tech- 
nique. 
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